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Ceramic Matrix Composites (CMCs) are 
Advantageous for High-Temperature Applications 



Typical higher operating temperature 
limits as compared to advanced 
metals 

Lightweight, tailorable properties 

Potential candidate materials in many 
aerospace structural applications 
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Objective of Present Work : Examine and Assess 
Approaches for Modeling Woven CMCs 



Material description - 

- Fiber is woven in a fabric preform of desired 
architecture such as plain weave, five- 
harness satin weave 

- Woven fabric is stacked in multiple layers 

- An interfacial coating is deposited on all fiber 
surfaces usually by a chemical vapor 
infiltration (CVI) process 

- The preform is then infiltrated by a silicon 
carbide matrix using the CVI process 

- The remaining porosity is filled either 
continuing to deposit the SiC matrix uisng the 
CVI process or in some cases using a slurry- 
cast melt-infiltration (Ml) process. 
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Multiple Scales in Woven/Braided Composite 
Continuum Modeling 
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Given Constituent Properties, Determine Effective Tow 
Properties 


Note : Porosity accounted for 
in SiC matrix properties 

Fully Dense E SiC * 420 GPa 


Correlated 

iBN Sylramic 

CVI-BN 

Intra-Tow 

Inter-Tow 

Constituent 

Fiber 

Coating 

CVI-SiC 

Ml-SiC 

Properties* 

0 

01 

c\i 

1204 °-C 

0 

01 

c\i 

1204 °-C 

0 

01 

c\i 

1204 °-C 

21 °-C 

1204 °-C 

Modulus (GPa) 

380 

365 


21 

14 


380 

358 

310 

276 

Poisson’s ratio 

0.17 

0.17 

0.22 

0.22 

0.17 

0.17 

0.17 

0.17 

CTE (10 6 / Q C) 

4.6 

8.0 

5.2 

10 

4.6 

9 

4.7 

9 

Th. Cond. (W/mK) 

43 

21 

3.1 

1 

70 

33 

68 

25 


*Ref: Murthy, P.L.N., Mital, S.K., and DiCarlo, J.A., NASA /TM- 1999-209 173 + in-house GRC properties 
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Analysis Methods Vary Based on Capabilities, Fidelity, 
and Computational Expense 



• Micromechanics Based (Analytical) Methods - 

- Laminate Approximation 

- W-CEMCAN - Generalized Method of Cells (GMC) 



CV1-S*C ♦ porosity • 


W-CEMCAN 


*Ref: Mital, S.K., Murthy, P.L.N., and Chamis, C.C., NASA/TM- 
1996-107165. 



Generalized Method of Cells 


1-Step and 2-Step Methods 


Ref: Aboudi, J.,Mech. Adv. Materl. Struct. 11, 329-366, 2004. 


• Also have High Fidelity Generalized 
Method of Cells (HFGMC) 

• Accounts for tow directionality 

• Elliptical tow shape 

• Account for tow undulation 
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Analysis Methods Vary Based on Capabilities, 
Fidelity, and Computational Expense 



• Hybrid Methods — 

User friendly 

Models many types of 2-D and 
3-D architectures 



pcGINA 

Ref: Gowayed, Y.A, Ph.D. Thesis, N.C. State University, 1992. 


• Numerical Methods (2-D and 3-D Finite Element Analvses)- 


Accurate 
Inefficient and 
expensive 



3D Finite Element Analysis (Abaqus) 
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Elastic Properties at Room Temperature 
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25 

CTE X (10' 6 / 9 C) 

TBD 

TBD 

4.63 

4.64 

4.66 

TBD 

TBD 
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*For 4 Load Cases 

Elastic properties within ~7% 

Constituent properties known no better 
Experimental repeatability no better 

Note : Execution times based on Intel dual core X7900 @ 2.8 GHz, 4 GB RAM 
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Efficient Methods Still Give Reasonably Accurate 
Approximations of Local Stress Fields 

• Von Mises stress fields predicted by models 

- 1204 Q C, applied in-plane strain of 0.1% 



Multiscale Lam. Theory, Rect. Tow - 9 2D subcells per ply - 0.015 seconds * Loading Direction 
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Multiscale Lam. Theory, Cross Tow - 25 2D subcells per ply - 0.015 seconds 


1-Step GMC-3D - 400 3D subcells - 0.08 seconds 


n 


HFGMC-3D - 400 3D subcells - 2.9 seconds 


HFGMC-2D - 656 2D subcells - 0.22 seconds 
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ImMAC Software Suite 

Integrated multiscale Micromechanics Analysis Code 



• Consists of MAC/GMC (stand-alone), FEAMAC(implemented within 
ABAQUS) and HyperMAC (implemented within HyperSizer) software 
codes. 


• ImMAC is software released by NASA GRC/RXL for multiscale 
analysis of composite structures 

• It links the behavior of a structure to the behavior of the composite 
constituent materials 


The key link between the scales is micromechanics , which provides 
the composite response based on the constituent behavior 

Constituent Composite Structural 

Response Micromechanics Model 

Repeating Unit Cell 


Fiber Interface Matrix 
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Creep Modeling 



• Multiscale Lamination Theory 

- Rectangular tow 

- Cross-shaped tow 



• 2D FEA 



- Transversely isotropic power law creep model used to model creep of tows 

- Five material parameters for transversely isotropic tows and three parameters for 
isotropic matrix 

- These parameters must be backed out from available measured composite creep 
curves 


Creep Parameters of Tows and Matrix Not Known 
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Efficiency is Critical when Backing Out Nonlinear 
Model Parameters 



• Typical to run hundreds of cases to determine parameters 

- Used efficient Multiscale Lam. Theory approach to obtain parameters 

• Utilized creep test results at 1315 Q C 

• Isotropic Ml SiC: k 7 = 55 MPa, ji = 3.7x1 0 10 MPa-s, n = 1 , rj = 1 , co = 1 

• Trans. Iso. Tow: k 7 = 6.9 MPa, ju = 6.9x1 0 10 MPa*s, n = 3.5, r| = 5, co = 5 


Experimental Creep Data, 1315 -C 



Characterized Tow and Ml SiC Creep Response 



Creep Time (hours) 
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Multiscale Lam. Theory Captures Primary and Secondary Creep, 
Even Though Constituents have Only Secondary 



• Stress redistribution from relaxation drives apparent primary creep 

• Effect of tow shape representation much more pronounced 

- Recall E affected only by ~0.5 % 
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Methods Give Similar Macroscopic Creep 

Predictions 
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Multiscale Generalized Method of Cells (MSGMC) For 
Concurrent Analysis of Woven/Braided Composites 
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Multiscale Generalized Method of Cells(MSGMC) 

Overview 



• Newly developed recursive GMC methodology 

• Each length scale in each subcell can call a separate GMC analysis 


• Works for any arbitrary multiphase material 
• Elastic / Inelastic / Damaae 
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Study Effects Of Micro, Meso, And Macro Parameters 

on Macroscale Response 
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Architectural Parameter 

Relevant Length Scale 

Values 

Tow Fiber Volume Fraction 

Meso 

0.46,0.48,0.50 

Tow Void Volume Fraction 

Meso 

0.01,0.05,0.07 

Tow Aspect Ratio 

Macro 

8,10,12 
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Range Of Macro Response Curves Given the 27 
Variations In Architectural Parameters 

Utilized Localized Void Model 



Architectural Variations 
clearly contribute to 
variation in measured 
material response. 

• Initial Modulus » 24% 

• UTS *2% 

• 1 st matrix cracking » 
16% 

• Post matrix cracking 
Modulus » 24% 

• 8 f impacted » 16% 
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Conclusions 


• Elastic Response - 

- Various methods for modeling 5 harness satin weave SiC/SiC 
composite were examined. 

- Methods generally fall in three categories: analytical, hybrid and 
numerical. 

- All methods do reasonably good job of predicting elastic 
properties as well as elastic stress fields. 

- Computational efficiency is the discriminator as analytical 
methods are orders of magnitude more efficient than fully 
numerical methods. It is important for multiscale analysis. 

• Creep Behavior - 

- Two methods (multiscale laminate analyses and 2-D FEA) were 
examined for modeling creep behavior. 

- Efficient methods are needed for backing out creep model 
parameters. 

- Creep/relaxation process among long, tows, trans. tows, and 
matrix is complex and drives the response. 

- Both methods agree reasonably well with experimental data. 
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Conclusions (contd.) 

• Multiscale Analysis of Woven Composites - 

- Demonstrated that a synergistic analysis using the multiscale 
generalized method of cells (MSGMC) can accurately represent 
woven CMC tensile behavior (loading/unloading). 

- Failure mechanisms are captured via local continuum damage 
model. 

- Variations in Weave Parameters (micro, meso, and macro) appear 
to contribute to variation in measured material macrolevel 
response. 
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